o 
o 

(N 
> 

o 



Oh- 
I 

o- 



> 

m 
(N 



o 



X 



Determination of the upper and lower bound of masslimit of degenerate fermionic 

dark matter objects 

Gabor KupQ 

Weizmann Institute of Science, PO Box 26, Rehovot 76100, Israel 

(Dated: February 2, 2008) 

We give a gravitational upper limit for the mass of static degenerate fermionic dark matter objects. 
The treatment we use includes fully relativistic equations for describing the static solutions of these 
objects. We study the influence of the annihilation of the particles on this mass limit. We give the 
change of its value over the age of the Universe with annihilation cross sections relevant for such 
fermions constituting the dark matter. Our calculations take into account the possibility of Dirac 
as well Majorana spinors. 
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I. INTRODUCTION 

Existence of a considerable amount of cold dark mat- 
ter in the present Universe has been regarded as a robust 
ingredient in recent astrophysics and cosmology. Several 
hypothetical particles were suggested in order to solve 
the problem of the cold dark matter. Some of them are 
bosons: axions, Q-balls. . . The others are fermions: neu- 
tralinos, axinos, gravitinos, quintessino. . . (One can read 
a good summary about the non-barionic dark matter can- 
didates in [1].) 

We study some general features of the degenerate 
fermionic matter here. Degeneracy means that every en- 
ergy states are occupied, that is such a system is spherical 
symmetric and does not rotate. Our aim is to determine 
the maximal possible mass of static degenerate fermion 
balls. 

First we make a naive non-relativistic estimation be- 
cause the non-relativistic system has been studied in 
the literature for decades. Several papers were devoted 
to the non-relativistic degenerate fermionic dark mat- 
ter with or wihtout barionic component (see for example 
P, 3, 4, 5, 6, 7]). Papers whose topic is some aspect of the 
relativistic behaviour of these systems are rare ([3]). We 
devote the second part of this paper to giving an exact 
relativistic mass limit for these systems. In the third part 
of the paper we study the influence of the annihilation of 
the particles on this mass limit. 

In non-relativistic case we use polytropic equation of 
state, in relativistic case we do not use such approxima- 
tion but the exact equation which one can get for degen- 
erate fermions. 

If we assume that the dark matter consists of neutral, 
weakly interacting fermions, their spatial distribution is 
driven by gravity and Fermi distribution. kT <^ rriQC^ is 
valid for such kind of particles today, that is if the den- 
sity is sufficiently high, the degenerate (T = K) Fermi 
distribution approximates the real distribution well, k is 



the Boltzmann constant, T is their temperature, toq is 
their rest mass and c is the speed of light. One can write 
that 



ZTTpp^ /3h'^ = g/niQ, 



(1) 



where g is the mass density of particles, pp is their Fermi 
momentum and h is the Planck constant. If the particles 
are Majorana spinors, the z factor is 8. If Dirac spinors, 
then 16. 



II. NON-RELATIVISTIC CASE 

We dispense with barionic matter and deal only with 
one component dark matter. This assumption is not good 
for cosmological scales but probably suitable for smaller 
scales. We derive the equations of the system first. 
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We can obtain pp from ^: pf — h i^-;^) ■ The 

energy density for non-relativistic case is 



e{p)p^dp 



we can integrate it easily 
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|e is valid for ideal gases. We know the equation of 
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Navier-Stokes equations are very simple for spheri- 
cal symmetric matter distribution in hydrostatic equi- 
librium: 



dP 
dr 



dU 
dr 



9, 



(4) 
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where r is the distance from the origin and U is the grav- 
itational potential. 



One needs also the Poisson equation V^t/ = AttGq, 
where G is the gravitational constant. It can be reformu- 
lated for spherical symmetric case: 



£U 2dU , ^ 
-—- + --—= 47rGe. 
ar'' r ar 



(5) 



If one substitutes P in ([4]) from ([3]) and integ- 
rates it, then gets relation between g and U: g = 

[■^{Uq — U)] , where Uq is the constant of integration. 

With this the Navier-Stokes equation is: — ^^"^ ' + 

2 d{Uo 



-jii = -^^G {^f {Uo - ur/\ 



d^ - -"-V5K/ v-u -; • If wc intro- 
duce the following variables x = y^r (where a = 
^^n'^Gh-^m^) and tJj = x{Ua - U), then it becomes 
very simple: 






f 
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(6) 



It is the well known Emden equation for n — 3. When one 
solves this equation, has only one free parameter, --p-\x=o, 
because ip(fi) rnust be in order to avoid singularity at 
x = 0. 



notation for the integral and A for the constant factor 
excluded the mo- Then we have M — -^k(xp). We 

know the general scaling feature of this equation M ex ^ 

(see for example), that is h{xf) = -§-, where the 

numerical value of B is approximately 132.384. If the 
mass reaches the Schwarzschild mass limit, we get a black 
hole with Msch ^ A^ ^ ^ ^ ^^. After some 



2GVa' 



algebra one can write relation between Msch and ttiq 
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(7) 



If Msch is measured in solar masses and toq in GeV, it 
means Msch ~ 2.11^ in Dirac case and Msch ~ 2.98^ 
in Majorana case. 

IV. RELATIVISTIC CASE 

For spherical symmetric case the metric can be written 
in the form 



ds^ 



e^c^dt^ ~r^{de^ +sin' 



')-e^dr\ 



III. MASS LIMIT IN NON-RELATIVISTIC 
CASE 

In the non-relativistic case we can make a raw esti- 
mation of upper limit of the mass of such a fermion ball. 
The mass density can be calculated from the previous re- 
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lations g — j^ ( ^ j . If we put this into the common 
expression of total mass, we get the following 

M ^An g{r')r'^dr' => 

^0 



M 



S^TTGy^ml Jo 
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where rp and xp denote the radius of this object, where 
the mass density drops to zero. We can introduce k{xp) 



where v and A are the function of r and t. We study the 
static case here, that is i^ := ^{^)j ^ = M''') J^]- 

One can obtain the Einstein equations from this and 
the stress-energy tensor of the ideal gas for static case: 
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where P is the pressure and e — gc^ + 3P/2 is the energy 
density of ideal gas. Prime denotes the differentiation 
along r. One can derive an explicite ordinary differential 
equation for i/" from these equations. 
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V = ::^z^ (8) 
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The mass density is the following: 



{mlc^pp + 2p%) ^Jml+p%/c^ - m'^c^l 
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The pressure can be obtained by means of integrating of the current of hnear momentum: 
I 
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We can not derive an equation of state with simple 
(polytrop) form from ([9]) and pO)) as in non-relativistic 
or highly relativistic limit. The best choice is to use pp as 
an independent variable instead of P or g. One needs an 
equation for pp now. It can be derived from the feature 
of stress-energy tensor T^^. = 0: 



dg 



(11) 



If one substitutes P and g in (fTTj) . he obtains for pp: 



p'f 



3h^ gippy + %P{pp] 
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ml+pl/c^iy' (12) 



The system ([8]), ([12]) contains two unknowns and two 
ordinary differential equations, that is it can be solved. 
((HI), (jlOp are only implicite auxiliary equations in tech- 
nical sense.) The calculations themselves follow closely 
Oppenheimer and Volkov's calculations [11|. 

This system has a scaling feature which can be seen 
without solving the equations. If gi and iy[ are solution 
of the system with r/ii, then the following functions also 
solve the system with 7712: 



92 = gi 



with the transformation of ri 



r2 = n 



m2 
mi 
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mi 
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A consequence of this feature is the following: 

M2 = Mi{^'' 
\m2 



(13) 



We see that the total mass is in inverse proportion to 
the square of the rest mass of the particles even as in 
non-relativistic case. 

In order to solve these equations we need initial condi- 
tions, v does not occur explicitely in the equations, only 
its derivatives. (Non-static equations contain v itself.) It 
means that its initial value is arbitrary. v'(0) is 0. This 
choice ensures that there is no singularity in the origin. 
One gets the value of p^(0) from t^'(O) by means (J12p . 
Only the central Fermi momentum, pp{0), remains free 
parameter. 

We can see the M{r) function of solutions for differ- 
ent Pp in FIG[TJ Solid curves belong to the relativistic 
case and dotted curves to the non-relativistic one. (We 
do not depict relativistic and non-relativistic curves with 
same central pp for sake of better perspicuity but we 
have tested that relativistic curves converge to the non- 
relativistic ones a pp{0) — *■ 0.) 

As one increases the value of pp{0), the end point of 
the non-relativistic curves start from the right bottom 
corner of the FIG[l]and evolve to the left upper corner. 

FIGlD shows the envelope of the end points of the 
curves in the FIG[TJ In non-relativistic case one can get 
solutions for arbitrary high masses but not in relativistic 
case. We can observe a "turning point" in the relativistic 
curve. There is a maximum mass here above which there 
are not solutions of the system. We denote this maxi- 
mum mass with M^ax and its location with Rmax- The 
ratio Mmax/Msch is ~ 5.29. 

(One can observe an interesting spiral structure in this 




FIG. 1; Comparison of relativistic and non-relativistic calcu- 
lations for mass as function of the distance from the centre 
for different central Fermi momentum, mo = lO^^eV, Dirac 
case. 



plot. We followed this "whirl" till the numerical resolu- 
tion allowed. This structure means several solutions with 
same mass. May be infinite for one M value if the spiral 
makes infinite number of revolution. Of course we can 
state similar for R.) 

FIG [3] shows the difference between the density distri- 
bution of such systems in Dirac and Majorana case. We 
chose one of the Dirac curves ending in a given point 
f^j^Dtrac^j^Diracy rpj^^j.^ jg ^^ Majoraua curvc which 

ends in same point but we could depict two corresponding 
Majorana curves. One of them ends at i?^*''''^ and the 
other at M^^^"-"^. The curve ending at M^*'""'^ spreads to 
larger volume than its Dirac counterpart because there 
is no so many energy states in the phasespace as in Dirac 
case, that is the central density is lower in this case. The 
curve ending at R^'^'^"-'^ contains more mass than its Dirac 
counterpart because this higer mass can " press" the Ma- 
jorana matter to the smaller volume. It is the qualitative 
reason why the Majorana curve in FIGj^] is shifted to 
the right and up. Quantitatively one can get the Majo- 
rana curve in FIG [2] from the Dirac one if r and M are 

multiplied by a factor ^J ZOirac/ ZMa^orana = \/2. 

If a fermionic ball reaches the Rmax with M > Mmax, 
then there is no static soultion, that is the gravitational 
collapse is unavoidable. We can conclude that this M„iax 
is the mass limit which we have been looking for. If we 
take into account the scaling feature, we have Mmax ~ 





FIG. 

mass. 



2: Radius of surface of our objects against their total 
mo = 10"eV. 



FIG. 3: Comparison of the density profile of Dirac and Ma- 
jorana case, mo — lO^'^eV. 



0.39^ in Dirac case and Mmax ~ 0.55^ in Majorana 
case (M is in units of solar mass and toq in GeV). 



V. INFLUENCE OF THE ANNIHILATION ON 

THE MASSLIMIT - LOWER BOUND OF THE 

POSSIBLE MAXIMAL MASS 



There are several papers about the annihilation cross 
section of weakly interacting particles which can be can- 
didates of dark matter. One can read a good review 
about this topic in Bertone's report 12] (see also the ref- 
erences therein) . We study the effect of the annihilation 
on the upper mass limit now. 

We can expand the velocity averaged cross section in 
a Taylor series: (av) = a + hv^ + . . . We do not know 
the value of the factors. One can only estimate upper 
limits of the total {av) in a given range of toq. We use 
Beacom's estimations [ij]. We calculate the final mass of 
such objects after 13.7 Gyear (age of the universe) with 
these upper limits started from the Mmax obtained in 
the previous section. It gives lower bound on the limit 
mass. (We have to note that it is not a lower bound on 
the mass of these objects because it is, of course, zero.) 

We do not make big error in the calculations of the de- 
caying time and the final mass probably if we use the 
{av) = constant non-relatvistic approximation. It is 
known that M(t)/Ma = . ^ in non-relativistic case 

[7|, where Mq is the initial mass and t is the character- 
istic time of the decaying (r oc Mq~^). It means that the 
decaying is very fast at the beginning and slower later. If 
we used higher terms in the expression of cross section, 
we just would accelerate the first phase but it would not 
influence the speed of the latter, non-relatvistic decay 
phase. FIGlDshows the ratio of the maximal velocity ob- 
tained from pp to speed of light against the ratio of the 
mass to the maximal mass from relativistic calculation. 
One can see that the maximal velocities of the particles 
decrease fast with decreasing M . We will see that the fi- 
nal mass is far below the hundredth of the original mass. 
It justifies our approximation. 

One can write the mass loss by the annihilation with 




FIG. 4: Ratio of the maximal velocity calculated from pp to 
speed of light against the ratio of M to Mmax ■ 




FIG. 5: Possible maximal mass M of the Fermion balls against 
the rest mass mo of Fermions. 



these assumptions in the following form: 



dM 
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32n(av) 
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gnr^dr = -{av)BM^ 
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This equation can be integrated easily: 

Mo 



M{t) = 



v/1 + 2{av)M^Bt 



(14) 



(15) 



If we use the scaling feature P^ . start the mass from 
Mmax and use Beacom's upper limits, then we get the 
following: 



M* 



M{t,mo) = 



\"io J 



1 + 2{(Jv)upper{mo)M,^ 



Bt 



(16) 

where M^^^ is the maximal mass at a given ttiq and 
{(Tv) upper (jno) dcnotcs that {av) depends on ni. If we set 



t = 1.37- 10"'^° year, then we can get the recent lower limit 
on the maximal mass of these Fermion balls. 

We summarize our result in FIGlH The upper two 
curves show the upper limit of the limit mass if there 
is not annihilation. The lowest two curved lines show 
the mass limit after 1.37 • 10^° year annihilation started 
from the upper curves, two curved lines above these two 
belong to 1.37 • 10® year annihilation. We depict them 
because Beacom's estimations of {av) give values for the 
cross section in the late universe. One can observe that 
even after this shorter time the mass limit is well below 
the limit without annihilation. We know further that if 
the {av) is larger than the natural limit, two horizontal 
lines, then the dark matter can not be thermal relic. 

Unfortunately the lower limit not strict but we can 
state this now. Probably the real cross section is much 
below the limits given by Beacom. If someone can give 
more strict upper bound on the cross section, our lower 
limit would increase. 



VI. CONCLUSIONS 

We studied self-gravitating systems of weakly inter- 
acting fermions. We determined the limit mass above 
which systems undergo gravitational collapse. We could 
see that this limit mass is a function of the rest mass of 
the particles: M = constant / m? . However, this relation 
was known before we give the exact value of this constant 
which depends on the type of fermions, namely whether 
they are Dirac or Majorana spinors. In Majorana case 
the central density cannot be as high as in Dirac case at 
same total mass. One can conclude from this that the 
value of this constant is higher in Majorana case than in 
Dirac case. It is an upper or absolute masslimit of such 
objects. 

We calculated the lower limit of the limit mass used 
the recent upper bound on annihilation cross section of 
dark matter particles. We could see that this limit is not 
strict but we can not state more on the base of the recent 
knowledge. 

Further investigations are needed to answer the ques- 
tion what is the probability that such objects gather mass 
above the mass limit and undergo gravitational collapse 
or how much dark matter they could " burn" if they could 
not collapse. 
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